Theoretical Study of Half-Doped Models for Manganites: 
Fragility of the CE Phase with Disorder, Two Types of Colossal Magnetoresistances, 
and Charge- Ordered States for Electron-Doped Materials 



H. Aliaga, 1 D. Magnoux, 2 A. Moreo, 1 D. Poilblanc, 2 S. Yunoki, 3 and E. Dagotto 1 

'National High Magnetic Field Lab and Department of Physics, Florida State University, Tallahassee, FL 32310. 
2 Laboratoire de Physique Theorique CNRS-FRE2603, Universite Paul Sabatier, F-31062 Toulouse, France. 
3 International School for Advanced Studies (SISSA), via Beirut 4, 34014 Trieste, Italy. 

(Dated: February 6, 2008) 

A comprehensive analysis of half-doped manganites is presented using Monte Carlo simulations applied to 
the double-exchange model with cooperative Jahn-Teller lattice distortions in two dimensions. A variety of 
novel results are reported. In particular: (i) The phase diagram is established in the A-Jaf plane, with A the 
electron-phonon coupling and Jaf the antiferromagnetic exchange between classical t2 g spins. The results 
include standard phases, such as the CE-insulating and FM-metallic regimes, but they also include novel states, 
such as a ferromagnetic charge-ordered (CO) orbital-ordered phase originally predicted by Hotta et ah. This 
state is compatible with recent experimental results by Loudon et al. ( ii) For realistic couplings, it was observed 
that the charge disproportionation S of the CO phase is far from the widely accepted extreme limit 5=0.5 of a 
3+/4+ charge separation. A far smaller 5 appears more realistic, in agreement with recent experiments by Garcia 
et al. and Daoud-Aladine et al. (Hi) Colossal magnetoresistance (CMR) effects are found in calculations of 
cluster resistances using the Landauer formalism. This occurs near the ubiquitous first-order phase transitions 
between the insulating and metallic states. The present result reinforces the previous conjecture that CMR 
phenomenology exists in two forms: the low-temperature CMR addressed here and the more standard CMR 
above the Curie temperature, (iv) The CE-state is found to be very sensitive to disorder since its long-range 
order rapidly disappears when quenched-disorder is introduced, contrary to the FM state which is more robust. 
This is also in qualitative agreement with recent experiments by Akahoshi et al. and Nakajima et al. (v) The 
phase diagram in the half-doped electron doping regime is briefly discussed as well. A charge-ordered state is 
found which is the analog of the x=0.5 CE phase. It contains a 3+/2+ charge arrangement at large A. Numerical 
results suggest that an approximate symmetry exists between the hole- and electron-doped systems in the large 
Hund coupling limit. 

PACS numbers: 75.50.Pp,75. 10.Lp,75.30.Hx 



I. INTRODUCTION 



Manganites are currently attracting considerable attention 
mainly due to the presence of the colossal magnetoresistance 
effect in magnetotransport measurements ll2 . In addition, these 
materials have a complex phase diagram with a plethora of 
ordered phases, a typical characteristic of correlated elec- 
tron systems. Many experimental and theoretical investiga- 
tions have unveiled the inhomogeneous character of the states 
of relevance to explain the CMR phenomenon, with a com- 
petition between ferromagnetic and antiferromagnetic states 
that induces coexistence of clusters, typically with nanometer- 
scale sizes^. The rationalization of this phenomenon, and 
concomitant explanation of the CMR effect, is based on the 
first-order transitions that separate the metallic and insulat- 
ing phases in the clean limit (i.e. without disorder)!^. The 
first-order character of the transition is caused by the differ- 
ent magnetic and charge orders of the competing states. The 
clean-limit phase diagram is illustrated in Fig. ^(upper panel). 
When quenched disorder is introduced in the coupling or den- 
sity that is modified to change from one phase to the other, 
the temperature where the Neel and Curie temperatures meet 
is reduced in value and eventually reaches zero as in Fig. ^ 
(middle panel). Upon further increase of the disorder strength, 
a spin disordered region appears at low temperatures, with a 
glassy behavior involving coexisting clusters (Fig. \l\ (lower 



panel)). Simulations by Burgy et al£ have shown that the 
clustered state between the Curie temperature and the clean- 
limit critical temperature (T*), with preformed ferromagnetic 
regions of random orientations, has a huge magnetoresistance 
since small fields can easily align the moments of the ferro- 
magnetic islands, leading to a percolative conductor in agree- 
ment with experiments 4 . The quenched disorder simply trig- 
gers the stabilization of the cluster formation, but phase com- 
petition is the main driving force of the mixed state. It has 
been speculated^ that similar phenomena should occur in a 
variety of materials, including, for example, the high temper- 
ature superconductors where nanoclusters have been found 5 . 



In spite of the phenomenological success of the previous 
calculations discussed above, there are many issues that must 
be further investigated in the manganite context to reach a rea- 
sonable understanding of these materials (for a recent discus- 
sion on open issues see Ref.|y). It is the purpose of this paper 
to contribute to manganite studies by addressing a variety of 
open topics, using computational techniques applied to real- 
istic models for Mn-oxides. In particular, here the phase dia- 
gram at hole-density x=0.5 is presented for the case of coop- 
erative Jahn-Teller phonons in interaction with carriers. Pre- 
vious studies mainly focused on non-cooperative phonons-. 
The phase diagram reported below in Fig. [2] has interesting 
features. For example, it contains the realistic FM charge- 
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FIG. 1 : Top General phase diagram of two competing phases in the 
absence of quenched-disorder (or when this disorder is very weak). 
Thick (thin) lines denote first (second) order transitions. Shown is a 
tricritical case, but it could be bicritical or tetracritical as well, g is 
some parameter needed to change from one phase to the other. Mid- 
dle With increasing disorder, the temperature range with first-order 
transitions separating the ordered states is reduced, and eventually 
for & fine-tuned value of the disorder the resulting phase diagram 
contains a quantum critical point. In this context, this should be a 
rare occurrence. Bottom In the limit of substantial quenched disor- 
der, a window without any long-range order opens at low temperature 
between the ordered phases. This disordered state has glassy charac- 
teristics and it is composed of coexisting clusters of both phases. The 
size of the coexisting islands can be regulated by the disorder strength 
and range, and by the prox imity to the original first-order transition. 
For more details see Refs. 20 The new scale T* discussed in the 
text- remnant of the clean-limit transition - is also shown. 



established the bicritical form of the £=0.5 phase diagram-. 
Moreover, a novel phase is found with both charge and fer- 
romagnetic order, in agreement with previous calculations^ 
and also with recent experiments by Loudon et alM Another 
feature of the results found here that is in good agreement 
with recent experiments (to be discussed below) is the value 
of the charge disproportionation. The widely accepted view of 
the extreme 3+/4+ separation needs to be revisited. A much 
milder charge separation appears more realistic. 

The x=0.5 phase diagram of the realistic model studied here 
is also analyzed in the presence of quenched-disorder, to test 
the results of phenomenological models, such as those shown 
in Fig. [2 Overall, good agreement with previous studies is 
found, but surprises are also reported such as the strong sen- 
sitivity of the CE phase to quenched disorder. This behav- 
ior is also in good agreement with recent experimental results 
by Akahoshi et al*^ and Nakajima et alM, where an asym- 
metry between the behavior of the FM- and CE-phases was 
observed when disordered was introduced in some materials 
with x=0.5. Finally, since electron doping of manganites has 
also been experimentally actively pursued 1 ', the model stud- 
ied here is analyzed for the case of electron doping as well, 
with an electronic density per site equal to 1.5. A state quite 
similar to the CE phase is observed, which could be found in 
future experiments. 

The overall conclusion is that theoretical studies of realis- 
tic models for manganites using unbiased techniques are un- 
veiling a remarkable qualitative agreement with experiments. 
While previous phenomenological descriptions of manganites 
are reasonable approximations for the understanding of the 
puzzling CMR phenomenon 2 , many other features observed 
in the present studies still need further investigations, such as 
the strong sensitivity to disorder of the CE state. 



II. METHOD AND DEFINITIONS 

The Hamiltonian studied in this paper is 

H = - ^ ^ 7 'K>r d i+a7'<T +h.C.) - JnJ2 Si ' Si 
ia77'cr i 

+ Jaf ^2 Si • Sj + A^(QiiPi + <52iT xi + Q 3 iT z i) 

+ (1/2) + + (i) 



disordered and CE charge-ordered phases well-known from 
experiments. In addition, further evidence is provided that 
the transition from the FM to the AF state is first-order at 
low temperatures. This first-order property allows us to ob- 
tain a huge magnetoresistance effect even working on small 
clusters, since small magnetic fields can unbalance the ground 
state from antiferro to ferro, in the vicinity of the competition 
region. For the calculation, techniques borrowed from meso- 
scopic calculations using the Landauer formalism are used. 
The shape of the phase diagram with increasing temperature 
is in good qualitative agreement with recent experiments that 



where di aa (diboO annihilates an e g -electron with spin er in the 
d x 2_ y 2 (d 3z 2_ r 2) orbital at site i, and a is the vector connect- 
ing nearest-neighbor (NN) sites. The first term is the NN hop- 
ping of e g electrons with amplitude t* , between 7- and 7'- 

orbitals along the a-direction: t* a =— V3t* h = — s/3t* =3t£ h =t 
for a=x, and t^ a =V3i^ h = V3t¥ =3t? h =t for a=y (our study 
will be restricted to two-dimensional lattices). Hereafter, t is 
taken as the energy unit. In the second term, the Hund cou- 
pling J H (>0) links e g electrons with spin s i= Y. ia p Aw*^ 
rfi7/3 (cr=Pauli matrices) with the localized t2g-spin Si, as- 
sumed classical with I Si 1=1. Jh is here considered as infinite 
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or very large. The third term is the AFM coupling Jaf be- 
tween NN < 2g spins. The fourth term couples e g electrons and 
MnC>6 octahedra distortions, A is a dimensionless coupling 
constant, Qn is the breathing-mode distortion, Q21 and Q^i 
are, respectively, (x 2 ~ y 2 )- and (3z 2 — r 2 )-type JT-mode dis- 
tortions, Pi = £ 7 (T d\ 1CT di ia , r xi = Y,A d L* d ib* +d\ ha d ia(y ), 

and T zi = Eo-C^W^iaa -4b,r d ib^)- The fifth term is the 

usual quadratic potential for adiabatic distortions and (3 is the 
spring-constants ratio for breathing- and JT-modes. In actual 
manganites, j3~2 (see Ref. fl4l) . and this is the value, unless 
something different is stated, we will consider throughout this 
paper. In undoped manganites, all oxygens are shared by ad- 
jacent MnOg octahedra and the distortions are not indepen- 
dent, suggesting that the cooperative effects are very impor- 
tant. This observation is likely valid even at finite hole densi- 
ties, since experiments show the presence of orbital ordering 
in the half-doped regime. To consider this cooperation, here 
oxygen ion displacements, denoted by uf and u\, are directly 
optimized™. 

It has been clarified in previous literature (e.g., Ref.0) that 
a large Hund coupling as used here suppresses double occu- 
pancy of the same orbital, and in this respect behaves as a 
Hubbard U interaction. In addition, a robust coupling A also 
acts as a U repulsion between electrons in the same site but 
different orbitals. As a consequence, it is a reasonable approx- 
imation to neglect the Coulombic interactions in the problem, 
which simplifies enormously the computational effort. 

The model will be analyzed primarily using a classical 
Monte Carlo (MC) procedure for the localized spins and 
phonons, in conjunction with exact diagonalization of the con- 
duction electron system. This last part of the process cor- 
responds to the solution of the single-electron problem with 
hoppings determined by the localized spin configuration. The 
resulting energy levels are then filled with the number of elec- 
trons to be studied, namely the simulations are carried out in 
the canonical ensemble (see Appendix). However, simula- 
tions directly in the grand canonical ensemble varying chem- 
ical potentials were also carried out. In the present work, all 
the calculations were made on a 2D 4x4 cluster, mainly with 
filling x=l/2, as explained in the introduction. Currently, it is 
not possible to comprehensively study larger systems, unless 
in special cases, due to the considerable CPU time that the di- 
agonalization process needs, and the large number of param- 
eters that must be varied to fully explore the phase diagram. 
In spite of the size limitation of our effort, the results reported 
here have a clear physical interpretation and size effects ap- 
pear to be mild for the quantities that were investigated. 

Updates of the spin and phononic 4>- u uf, uT\ config- 
urations are accepted or rejected according to the Metropolis 
algorithm. The number of MC steps per site is typically taken 
as 3000 for thermalization, with an additional 10000 for mea- 
surements. The simulations usually started with random states 
at high temperature, and then the temperature was decreased 
slowly, but at very low temperatures some of the simulations 
started with ordered states to speed up the convergence. The 
results of the simulations were analyzed using three separate 



but related quantities: the spin structure factor 

5(k)=^(S i .Sj)e ik ( r - r j), (2) 

of the classical spins, the charge structure factor 

iV(k)=^(p i p j )e ik -( r - r j ), (3) 

ij 

and the orbital structure factor defined as, 

T(k) = ^(Ti.Tj)e ik -( r '- r j), (4) 

ij 

with T ; = (r xi ,r y i,r zi ) and r yi = i Y,a( d L<r d ibc -d\ ha d- laiJ ) . 
Frequently, visual investigations of snapshots of the spin and 
phononic configurations at low temperatures were useful to 
guide the intuition. Some dynamical properties were also 
studied, such as the density-of-states. For details see Ref. 0. 
More sophisticated techniques used to study the conductance 
of clusters are described later in the paper. 

HI. LOW-TEMPERATURE PHASE DIAGRAM 

The phase diagram of the double-exchange model with co- 
operative JT phonons at a;=0.5 is shown in Fig. [2] The proper- 
ties of each phase were analyzed through the spin, charge, and 
orbital correlations in real and momentum space defined in the 
previous section. Varying Jaf three spin regimes were iden- 
tified, similar to those reported in previous investigations 16 . 
Typical results for the evolution of the energy with Jaf are 
presented in Fig- clearly showing the three states that com- 
pete at fixed A. Using this procedure, the phase diagram was 
constructed. Note the clear first-order transitions separating 
the different phases, shown in the figure as level crossings. 

In agreement with well-known experimental results for 
half-doped real manganites, a CE phase is stabilized in the 
intermediate Jaf range. The actual correlations for this phase 
will not be shown explicitly here since they have appeared in 
some previous investigations^ and our results agree very well 
with those results. However, there are issues to remark not dis- 
cussed before. For example, note that the CE phase is present 
not only at large A but also at small electron-phonon coupling. 
This last result suggests that a nonzero Jaf is sufficient to sta- 
bilize the zigzag chains characteristic of the CE region. Hotta 
et aZJi proposed that a band-insulator picture is needed to un- 
derstand this phenomenon, and the present results support this 
view. It is the geometry of the zigzag chains that plays the key 
role in stabilizing the CE state. This result may have important 
consequences for experiments as discussed below. The dotted 
line represents a continuous transition between the standard 
CE phase and a novel canted CE' phase. To understand the 
characteristics of the new CE' state, note that the standard CE 
phase can also be interpreted as an arrangement of parallel AF 
zig-zag chains, displaced one lattice parameter in the x- and 
y-directions simultaneously. When the relative spin angle be- 
tween the neighboring AF zig-zags is 71", we have the normal 



4 



2.4 
2.2 
2.0 
1.8 
1.6 
1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



- 1 

- \ 


Y 


i i 


i 1 i 1 

2D ■ 


.FfX 








CO 








"00 






AF " 








OD 


- 


\ 










CE \^ 






FM \ 


CO \ 






i OO 






CD 








OFM 




CE' ^S." 






CO 








OO 

i 1 i 1 



0.0 



0.1 



0.2 0.3 0.4 0.5 



0.6 



FIG. 2: Phase diagram of the half-doped 2D two-orbitals model with 
cooperative Jahn-Teller phonons at T/t^0. The Hund coupling Jh 
is infinite. The phases were obtained analyzing the crossing of en- 
ergies (as explained in Fig. [5), and also S(k), T(k), and -/V(k) (see 
definitions in text). The magnetic phases present are FM, CE, CE' 
and AF. The symbols CO(CD) stand for charge-ordered (-disordered) 
phases. OFM corresponds to a homogeneous state with dominant 
uniform x 2 — y 2 orbital order. The staggered orbital-ordered states, 
OO, have the standard "CE" pattern, present in the three phases. The 
FM-CD-OFM phase is metallic, the rest of the phases are insulating. 
Thick lines denote first-order transitions, the dash-line is second or- 
der. The lines are smooth interpolations between a finite but large 
number of points obtained numerically. 



CE phase. In the CE' regime this relative angle varies contin- 
uosly as A or Jaf are increased. Finally, when this spin angle 
between neighboring AF zig-zag chains is 0, we recover the 
normal AF phase. 

At sufficiently large Jaf coupling, an antiferromagnetic 
spin arrangement is found. This phase has not been observed 
yet experimentally, but there is no reason to believe that the 
state cannot exist in some Mn oxide. Within the accuracy of 
our study, this AF phase is believed to have no long-range 
order in the orbital and charge degrees of freedom. This cor- 
responds to a nearly uniform arrangement of charge at small 
A. However, at large electron-phonon coupling the situa- 
tion changes qualitatively. In this regime, a "polaronic" state 
seems to be stabilized, where the 8 electrons present in the 16- 
sites cluster studied here become trapped in random locations. 
If in this random distribution some of the electrons are close to 
one another, local orbital order similar to that of the x=Q limit 
was observed. Again, this AF phase has not been observed 
experimentally yet, but the recent progress in the discovery 
of new phases (see next section) suggests that this remains a 
serious possibility. 

The regime of small Jaf coupling is dominated by FM 
states. The large-A zone has a novel state which is dis- 
cussed in the next section. At intermediate and small cou- 
plings, the state that dominates at low temperature has the 
characteristics of the FM metallic state well-known to exist 
in three-dimensional Mn oxides. For instance, in the charge 
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FIG. 3: MC energy vs. Jaf for x=l/2, A=1.8, T/t^0. This low- 
T energy is calculated after starting with a random high-T state and 
lowering T slowly. For this value of A, the three FM, CE and AF 



phases have OO with the d- c 
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usual alter- 
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nation of orbitals, distributed in zig-zags. The x 2 — y 2 orbitals at the 
"Mn 4+ " sites are populated at small and moderate values of A. The 
three phases also have checkerboard-type charge ordering. Near the 
crossings of levels we have observed a notorious slowing down in 
the convergence of the MC calculations. Lines are guides to the eye. 
Very similar results were obtained at several values of A and were 
used to construct Fig. [2] 



sector there is no order. Regarding orbital order, the corre- 
lations suggest the presence of a uniform arrangement. This 
uniform orbital order at small A is not produced through a 
spontaneous-symmetry -breaking process, but it exists as a di- 
rect consequence of the form of the kinetic energy in our two- 
dimensional simulations. To show this, the procedure is the 
following. Consider the limit of Hund coupling infinite as in 
our studies, and also A=Jaf=0. In this case, the Hamilto- 
nian can be diagonalized exactly in momentum space and the 
ground state for the 4x4 cluster can be exactly constructed. In 
this state, the mean value of the number operator for the two 
orbitals of relevance can be calculated as well. The results 
are not equal for the two orbitals, but there is an asymmetry 
in favor of the x 2 -y 2 orbital. Then, in this respect the uni- 
form orbital-order in the state is explicit in the model, and 
it is not induced by the JT phonons. These conclusion were 
also checked on larger lattices, such as 20x20, and in three- 
dimensional clusters. 



IV. NOVEL FM/CO PHASE 

It is interesting to remark that the phase diagram in Fig. [2 
contains not only the FM-metallic and CE-insulating phases 
- well-established experimentally - but other phases as well. 
For example, at large Jaf an AF phase exists in a wide range, 
and this phase could be observed in future experimental inves- 
tigations, as already discussed. Even more exciting is the case 
of the ferromagnetic charge-ordered (FM/CO) phase in the 



5 



upper-left corner of Fig. [2] This phase was found in previous 
investigations by Hotta et ali using cooperative phonons (and 
suggested by Yunoki et al& as well, using non-cooperative 
phonons). It is quite remarkable that recent experimental ef- 
forts by Loudon et ah 10 have unveiled the presence of a FM 
insulating phase at half doping, compatible with the FM/CO 
state found in the Monte Carlo simulations. For complete- 
ness, in Fig. |3 the orbital arrangement of the FM/CO phase 
found in our simulations is shown, although it has been al- 
ready presented in RefQ The orbital order corresponds to the 
same pattern of 3a; 2 — r 2 /3y 2 — r 2 orbitals of the CE-phase, 
but the spins are arranged ferromagnetically, since Jaf is not 
large enough to induce the CE order. Future experiments will 
clarify if indeed the theoretically predicted FM/CO regima^ 
corresponds to the results reported in Ref. Note that the 
discovery of more phases than previously known is a recent in- 
teresting trend, both in experiments and theory. In fact, at x=Q 
a novel E-phase has been observed in simulations by Hotta et 
as well as in experiments by Kimura et al*^. Many sur- 
prises may still be found in manganite investigations, even at 
the rather elementary level of characterization of states, since 
its strongly correlated character creates a plethora of phases 
that are in strong competition. 




FIG. 4: Orbital, spin, and charge arrangement for the FM-CO-OO 

o | 

phase. The 3x — r and 3y — r orbitals represent the Mn 
ions, while the empty sites represents the Mn 4+ sites. This state was 
presented for the first time in Ref.0- Evidence of a FM phase with 
checkerboard CO pattern has been found recentlyp- 



V. PHASE DIAGRAM VARYING TEMPERATURE 

The properties of the model studied here were also investi- 
gated at finite temperatures. In this case, the errors in the es- 
timations of the phase boundaries are larger than at zero tem- 
perature. The reason is that at T=0, the first-order transitions 
can clearly be established through level crossings, which exist 
even in small systems. However, the continuous transitions 
that are found at finite T can only be roughly located on small 
systems based on spin, charge, and orbital correlations at the 
largest distance available in the clusters studied. In spite of 



this problem, the results shown here are sufficiently accurate 
to understand the main trends in the phase diagrams. 

Typical results can be found in Fig. [5] The level-crossing 
procedure shown in (a) establishes easily the phase diagram 
at low T. At finite T, the momentum space correlations 
show a rapid increase upon cooling at characteristic temper- 
atures (see Figs. |5Jb,c,d) as typical examples) due to the in- 
crease of the real-space correlations at the largest available 
distances. These temperatures are the best approximations to 
the true bulk-limit critical temperatures that our investigations 
can produce at present. 
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FIG. 5: (a) MC energy vs. Jaf for 2=1/2, A=2.0, T/t^0. The loca- 
tion of the level crossings were used in the low-T phase diagram of 
Fig.|21 (b) Spin structure factor S(q) as a function of T/t, for q=(0,0) 
and Jap /t=0, indicating that at low T the FM phase is stabilized, (c) 
S(q) vs. T/t, for J A F/i=0.05, A=2.0, and momenta q=(0,7r) and 
q=(7r,0), characteristic of the CE phase, (d) Charge structure factor 
N(q) vs. T/t for the same parameters as in (c). In (b), (c), and (d) 
the approximate temperatures where correlations become robust are 
indicated by arrows. 

An example of the phase diagrams constructed by this pro- 
cedure is given in Fig.|6| where the results at A=2.0 are shown. 
Here, the 'oscillations' in the characteristic temperatures are 
indicative of the errors in our procedure, and of the clear ten- 
dency in simulations at large A to spend considerable Monte 
Carlo time trapped in competing states, as it occurs in glassy 
systems. In spite of these complications, the phase diagrams 
such as Fig.|6]are sufficiently informative to unveil the domi- 
nant properties of the system. 

Some of the properties of Fig. [5] are worth explicitly dis- 
cussing: (i) The Neel critical temperature of the CE state is 
the lowest among the three dominant states. This may ex- 
plain in part its sensitivity to disorder discussed later in this 
paper, (ii) At large A, charge order in the CE phase occurs at 
a temperature larger than the Neel temperature, as it occurs 
in many manganite experiments, (iii) In the range of Jaf be- 
tween 0.02 and 0.03, the CE phase is stable at low-T, but upon 
heating the FM state is stabilized, before the system becomes 
paramagnetic. This curious behavior is compatible with ex- 
perimental results recently reported by Tomioka and TokuraS, 
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FIG. 6: Phase diagram T/t vs. JAp/t of the 2D two-orbitals 
cooperative-phonon model at A=2.0. To, Tce, and Taf are repre- 
sented by squares, triangles, and asterisks, respectively. These tem- 
peratures are a crude estimation of the magnetic ordering tempera- 
tures of the FM, CE and AF phases, obtained in our small cluster sim- 
ulations, as explained in the text. The oscillations in the results are 
indicative of the error bars in the critical temperatures. The tempera- 
tures (Tco) where charge and orbital correlations become important 
upon cooling are marked by circles and a dotted line. Transitions to 
the paramagnetic (PM) and CO phases are second order. Transitions 
FM-CE, and CE-AF are first-order at low temperatures. 
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7: Experimental phase diagram of the half-doped compound 
o.sMnO-j,, reproduced from RefU The charge- 
and orbital-ordered insulator and ferromagnetic metallic states are 
denoted by CO-OOI and FM, respectively. The transitions from (to) 
the CO-OOI phase are represented by circles. The transitions to the 
FM phase are represented by squares. The Neel temperature is rep- 
resented by solid triangles. The grey area indicates hysteresis. Note 
the presence of a FM phase above a CO/AF phase, compatible with 
our results, as discussed in the text. 



reproduced in Fig. It is remarkable that the model stud- 
ied here is able to qualitatively reproduce even this fine detail 
of the real phase diagram of x=0.5 manganites. The differ- 
ent dimensionalities between experiments and the simulations 
reported here are not a problem, since it is well-known that 
the phase diagrams of the double-exchange model are quali- 
tatively similar in the dimensions of interest^. 

In Fig. [8] the numerically obtained phase diagram at A=1.3 
is also shown, as representative of the intermediate A range. 
As in the previous cases, the information was obtained on a 
16-sites lattice, using spin correlations at the largest available 
distances to decide which tendency dominates in the ground 
state. Instead of showing explicitly the numbers defining the 
lines as done in Fig. [6j here a smooth average is presented 
for simplicity. Error bars similar to those of Fig. [5] should be 
assumed at A=1.3 as well. In good agreement with the re- 
sults at other As, the Curie temperature in Fig. [8] is found to 
be larger than the CE critical temperature, and there is a re- 
gion of couplings where ferromagnetism is stabilized above 
the CE state, before turning paramagnetic upon further heat- 
ing. The gray region denotes a regime where charge-ordering 
appears, but not spin order (however, note that this region with 
phase competition is particularly complicated to study due to 
the presence of many competing minima in the free energy). 

Figure|9]illustrates the presence of gap and pseudogap (PG) 
features in the density-of-states (DOS), due to the robust value 
of A used. At low temperatures in the CE phase, there is 
clearly a gap due to charge ordering. However, the existence 
of a reduced DOS at the chemical potential in the form of a 
pseudogap survives the increase of the temperature, and above 
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FIG. 8: Phase diagram of the two-orbitals model with cooperative 
Jahn-Teller phonons at A=1.3 and /3=100. The results are smooth in- 
terpolations using the numerically available data, and in this respect 
they are not exact. The notation is standard. Note that at this cou- 
pling there are pseudogap features in the density-of-states above the 
ordering temperatures, and a hard gap below those critical tempera- 
tures. The gray area is a region where the simulations give a mixture 
of orders, for the large number of iterations carried out. 



the Neel temperature the effect appears to be caused by dy- 
namical Jahn-Teller distortions. This result reinforces the no- 
tion that a large depletion of the DOS at the chemical potential 
should be present in a broad range of couplings and tempera- 
tures in models for manganites. The existence of a PG feature 
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FIG. 9: Density-of-states obtained numerically at the couplings and 
temperatures indicated. Comparing with Fig. [8] these results show 
that pseudogap features appear in the DOS above the ordering Neel 
temperature of the CE state, at A=1.3 and ,5=100. Pseudoga ps have 
also been observed in many previous investigations (see Refs. l2fl2ll 
E3. for instance). 



severely questioned*! 

How can we understand this challenging new phenomenol- 
ogy? A recent key observation by Hotta et alJ-L indicates 
that the CE state may have an origin totally different from 
the "large Coulomb repulsion" view of Goodenough (or the 
large A analog in the presence of JT distortions). In the new 
perspective, the formation of the zigzag chains of the CE state 
can occur even at very small A since these chains lead to the 
optimization of the kinetic energy in the presence of a nonzero 
Jaf- In this context, the insulating character arises mainly 
from a band-insulator picture, already present in the study of 
individual zigzag chains. The validity of this view is clear in 
the phase diagram of Fig. |2] that shows a CE phase that ex- 
tends to A=0, a result quite difficult to understand from the 
large-A standard perspective. Since the FM metallic and CE 
phases are not in contact at large A in the theoretical phase 
diagram, and since experimentally they are known to be close 
in energy (at least for LCMO), then the 3+/4+ view of the CE 
state must indeed be revisited. 



was first noticed theoretically by Moreo et al. |20j and exper- 
imentally in photoemission experiments by Dessau et al*^. 
More recently, in optical conductivity experiments by Noh 
and collaborators^, PG features were also observed. Pseudo- 
gaps in DOS alter substantially the transport properties of the 
system, and in addition they suggest interesting similarities 
with other materials -such as the cuprates- where pseudogap 
features have also been identified. In fact, most of the phe- 
nomenology of the Mn- and Cu-oxides is similar, with phase 
competition being an important ingredient to understand phe- 
nomenologically their properties™. 



CO 




VI. CHARGE DISPROPORTIONATION 

Our calculations allow us to address theoretically a recent 
controversy in the experimental literature related with the ac- 
tual value of the charge disproportionation in charge-ordered 
half-doped manganites. The "standard folklore" since the 
1950s says that the CE state is made out of a checkerboard 
distribution of 3+ and 4+ charges. However, on theoretical 
grounds this state appears to be too "extreme". Since the 
work of Yunoki and collaborators 8 establishing the phase di- 
agram of manganite models with JT phonons at x=0.5, it has 
been noticed that the CE phase is close to a ferromagnetic 
metallic phase only at small and intermediate values of the 
electron-phonon coupling A. Having the two phases close to 
one another is important to understand the phenomenology of 
Lai-^Ca^MnOa (LCMO), where it is known that the two 
phases "touch" at a;=0.5, in the phase diagram varying the Ca 
concentration. If indeed A is not very large, then the charge- 
ordered state should be far from the 3+/4+ extreme separation. 
In fact, recent experimental work by Garcia et al, Daoud- 
Aladine et al, and others have revealed a very weak signal for 
charge disproportionation, and the standard 3+/4+ has been 



FIG. 10: Charge disproportionation S vs. A, calculated in the CE 
phase of Fig. [2] Starting at A=0, the values of Jaf used are 0.25, 
0.25, 0.25, 0.25, 0.2, 0.2, 0.2, 0.15, 0.15, 0.1, 0.1, 0.05,and 0.05, 
due to the tilted shape of the CE phase. A value S ~0.5 corresponds 
to the standard charge ordering, with Mn 3+ and Mn 4+ arranged in 
a checkerboard pattern. To reproduce this extreme charge-ordered 
state an abnormally large value of A is needed. More realistic param- 
eters correspond to a far less dramatic charge separation. 

Figure ^3 shows the numerically obtained charge dispro- 
portionation S along the CE phase varying A (note that the 
bended nature of the CE phase in the phase diagram Fig. [3] 
forces us to adjust Jaf as A changes). The entire CE region 
has a nonzero 5 according to Fig.^3- compatible with its 
insulating nature - but clearly there are three qualitatively dif- 
ferent regimes. At small A, the small value of 5 suggests that 
the charge ordering arises as a consequence of the more dom- 
inant zigzag chain formation, as predicted by Hotta et alll, 
namely the spin order dominates over the charge ordering. In- 
creasing A, a crossover regime is observed in Fig.^| which 
is itself followed by the extreme case of 3+/4+ separation at 
very large A. Since the phase diagram of Fig. [2] shows that 
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the FM charge-disordered and CE phases are in contact only 
below A=1.2, this implies an upper bound of approximately 
5=0.15 in the charge disproportionation, a result far from the 
more standard assumption (5=0.5. If A is slightly reduced from 
1.2 to 0.8, <5 can be as small as 0.05, compatible with the re- 
cent experimental information. As we have stated before, a 
large Hund coupling suppresses double ocupancy of the same 
orbital. This interaction together with the electron-phonon 
interaction and elastic energies behaves like the Hubbard U 
interaction™. The exact equivalence have not yet been clearly 
established. Our results are in good agreement with van den 
Brink et al. in the sense that turning on A in our study (U in 
Ref l25l) . the charge migrates from the corner sites of the zig- 
zag chain, to the bridge ones. However, the maximum charge 
disproportion available in our model (5=0.5) differs consider- 
ably from the (5=0.185 result obtained for U=oo. These dis- 
crepancies deserve further investigations to clarify the range 
of validity of both models. From all these considerations, our 
conclusion is that indeed the standard extreme view of the CE 
state needs revision. In fact, the surprising sensitivity of the 
CE state to the introduction of disorder (to be discussed be- 
low) may be indicative of a fragility which intuitively appears 
more related to the small-A zigzag-chain-driven regime than 
to the large A limit. Clearly, the last word on the nature of the 
CE-state has not been said, and more work should be devoted 
to this issue. 



VII. INFLUENCE OF QUENCHED DISORDER ON THE 
PHASE DIAGRAM AND FRAGILITY OF THE CE PHASE 



In Fig. II lh. the numerically obtained phase diagram at fixed 
A=l varying Jaf is shown in the absence of disorder (clean- 
limit). The shape is similar to results presented in previ- 
ous sections at other As, with competing FM and CE phases. 
Shown are also lines of constant spin-spin correlations, which 
will be useful to understand the results in the presence of 
quenched disorder. In Fig. II lb. the modified phase diagram 
after introducing disorder in Jaf is presented. The disorder 
is included by selecting the coupling Jaf between localized 
spins at every link as the uniform value plus an extra contri- 
bution AJaf, which is randomly added or subtracted (i.e. a 
bimodal distribution is used). In agreement with the expec- 
tations described in the introduction, Fig.Q including disor- 
der the lines of constant spin correlations bend and reach zero 
temperature at particular values of Jaf- This is to be con- 
trasted with the clean-limit case in Fig. II lh. where the lines of 
constant correlations 'collide' with the CE phase, rather than 
reaching zero temperature as with disorder incorporated. The 
behavior of Fig. II lb is indicative of the presence of a spin- 
disordered window separating the FM and CE phases, as in 

Fig.nt. 

However, the situation appears to be more complicated than 
previous discussions would have suggested. The reason is that 
the CE phase was here found to be abnormally sensitive to 
disorder. At least for the value of disorder used, a glassy state 
is obtained in approximately the same region where the CE 
phase was found to be stable in the clean limit. The glassy 
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FIG. 11: (a) Phase diagram T/t vs. Jaf / 1 for the clean limit, with 
A=l. The closed circles and diamonds indicate the transition tem- 
perature at which the long-distance correlation < SiSj V > approx- 
imately vanishes upon heating (see text). This is indicative of the 
temperature where long-range order develops. Dashed and dash- 
dotted lines represent the lines of constant < SiSj V >=0.2, and 
0.5, respectively, (b) Phase diagram T/t vs. Jaf /t for the disorder 
parameter AJaf=0.07, and A=l. The new symbols open -squares 
and -triangles represent the transition temperatures to a glassy state 
with short-range correlations but with no obvious long-range FM and 
CE order, respectively. The CE state was found to be very suscepti- 
ble to disorder, (c) < SiSj V > vs. Jaf ft for AJap=0 (squares), 
A Jaf=0.04 (triangles), and A Jaf=0.07 (circles) at low T. These re- 
sults are the average of five different configurations for the bimodal 
distribution J, = Jaf ± A Jaf in the AF links. 



nature of the CE region appears in S(k) and it is very clear 
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in Monte Carlo snapshots (not shown), that present distorted 
zigzag chains and patterns of charge ordering - nearly frozen 
as the Monte Carlo time evolves - in a glassy-looking arrange- 
ment. The effect is also manifested in the small value of corre- 
lations, as shown in Fig. II It which contains the spin correla- 
tions at the largest distance on the small cluster used. Without 
disorder, this correlation has the largest allowed value in both 
the FM and CE phases and as the disorder increases, the cor- 
relation reduces its value as expected. What is remarkable is 
that this effect occurs clearly more rapidly in the CE region 
than in the ferromagnetic phase. In fact, at A Jaf=0-07 the 
spin correlations are nearly negligible in the CE regime, indi- 
cating the disappearance of order at the largest distances here 
available. This is to be contrasted against the behavior of cor- 
relations in the FM phase, which are far more robust. 



The present results are in excellent agreement with the re- 
cent experimental studies of Ln^Lai^MnOs (Ln = rare- 
earth) by Akahoshi et a/jAI and Nakajima et al&. In these ex- 
periments it was observed that the phase diagram has bicritical 
behavior (Fig. II lh) when the crystals are grown slowly, mini- 
mizing by this procedure the influence of disorder. However, 
upon rapid quenching of the growing process to introduce dis- 
order explicitly, it was observed that while the FM phase re- 
duces its critical temperature by a reasonable amount, the CE 
phase is much more affected, turning into a glassy state at low 
temperatures. This experimentally-observed abnormal sen- 
sitivity of the CE phase to disorder -here nicely reproduced 
in the Monte Carlo simulations- is expected to be an impor- 
tant ingredient for the understanding of the phenomenology 
of manganites. Disorder indeed plays a key role in the phase- 
separated scenario for Mn oxides™, and here its relevance is 
further confirmed. The simple view expressed in Fig.^ with 
a symmetric behavior between the two competing phases, is 
more complex in practice. 



The disorder sensitivity of the CE-phase may arise from 
the fragility of the zigzag chains. As argued before, the realis- 
tic regime of couplings, with A of order 1, presents a charge- 
ordering pattern which is far from the extreme limit 3+/4+, 
rendering the staggered charge-ordering less robust than pre- 
viously believed. In addition, the one-dimensional nature of 
the zigzag chains also contributes to this fragility. Imperfec- 
tions can easily destroy these nontrivial geometrical arrange- 
ment of chains. Unfortunately, the limited size of the clus- 
ters used here does not allow us to investigate in more de- 
tail the observed CE fragility. These important issues should 
be analyzed in more detail in future theoretical and experi- 
mental investigations. It is interesting to remark that similar 
conclusions have been recently reached using a one-orbital 
model with cooperative phonons 26 . Reaching independently 
the same conclusion using different models, methods, and lat- 
tice sizes provides confirmation of the fragility of the charge- 
ordered state that competes with ferromagnetism. This is an 
important issue that deserves further studies. 



VIII. TWO TYPES OF CMR 

Let us start the analysis of d.c. transport properties in the 
half-doped system studied here, aiming toward a better un- 
derstanding of the CMR effect. In our opinion, these effects 
should actually be divided in two classes: (1) the CMR ef- 
fect that occurs at low temperatures and corresponds to very 
abrupt first-order transitions, and (2) the more standard CMR 
effect that appears in the regime immediately above the Curie 
temperature. Years ago, Tokura and collaborators already sug- 
gested these two types of CMRs 27 . Fi gurelT2"lprovides a direct 
experimental evidence of this phenomenon, with a huge effect 
at low temperatures and a relatively smaller effect near Tc ur i e . 
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FIG. 12: Temperature dependence of resistivity under various mag- 
netic fields for (Ndi_ y Smj / ) 1 /2Sr 1 / 2 Mn03 with y=0 (a) and 0.75 
(b). The hatched area represents thermal hysteresis. Results from 
Tokura et alw— 

This phenomenology - with two CMRs - has recently been 
discussed by one of the authors (E.D.) in Ref. 0, from where 
Fig-Elis reproduced. In this figure, the "CMR1" effect cor- 
responds to a direct field-induced transition from the AF to 
the FM phases, upon the application of an external field. In 
this case, quenched disorder is not needed and the transition 
occurs even in the clean limit. For CMR1 to occur, it is suf- 
ficient to be located on the insulating side but close to the 
first-order transition. Under these circumstances a relatively 
small field - which obviously favors the FM phase - renders 
ferromagnetism more favorable and an AF to FM first-order 
jump is induced. On the other hand, Fig. ^] also shows the 
standard "CMR2" proccess which is believed to originate in 
the more complex percolative regime induced by a nanoscale 
phase separated state above the Curie temperature*^. The re- 
alization that there are two types of CMRs substantially clari- 
fies the magnetotransport phenomenology of Mn oxides. 

To address explicitly the first type of MR transitions (i.e. 
CMR1), a clean-limit investigation should be sufficient. To 
study the resistance of a cluster vs. temperature in the 
presence of magnetic fields, a setup similar to those often 
used in investigations of mesoscopic systems will be here 
employed^. This method was recently applied to a techni- 
cally similar problem of spin and carriers in interaction in the 
context of diluted magnetic semiconductors 2 ** The interact- 
ing cluster under investigation is assumed connected to "ideal 
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FIG. 13: Schematic representation of the generic phase diagram in 
the presence of competing FM metal and CO/AF insulator, and for 
quenched disorder not sufficiently strong to destroy entirely the first- 
order transition at low temperatures, g is a generic variable needed 
to transfer the system from one phase to the other. CMR1 and CMR2 
are the regions with two types of large MR transitions, as described 
in the text (see also Ref.P 



leads", as shown in Fig.^J The information about these leads 
is included through exactly calculated self-energies which are 
located at the cluster boundaries. The current circulates after 
an infinitesimal voltage drop is included. The technical as- 
pects have been recently reviewed by Verges^ and they will 
not be repeated here. This technique appears to be better than 
other methods often used in numerical simulations, that rely 
on the analysis of the Drude weight D of the optical conduc- 
tivity with the cluster assumed in isolation. In this case, D 
is often negative due to finite size effects. In addition, the 
associated zero-frequency delta function must be given an ar- 
bitrary width to recover a finite resistivity. These problems are 
avoided in the present formulation. 



ideal 




ideal 


lead 


lead 



FIG. 14: Geometrical setup used here for the calcul ation of the resis- 
tance of a cluster. For more details see Refs, l28l29l 

A typical result obtained with the setup of Fig.^^is shown 
in Fig.^J To transform from resistance to resistivity, a cross 
section of size 4a x la (a=lattice spacing) is assumed. In the 
clean limit, and using as reference the phase diagram of Fig.|5J 
it was observed that for a range of Jaf the material behaves 
as an insulator since the CE phase dominates, while for other 
couplings a metallic behavior is found once in the FM regime 
(here insulator and metal are defined simply based on the sign 
of the slope of the resistance vs. temperature curves). The 
technique used here neatly reproduces the expected behavior 
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FIG. 15: Clean-limit investigation of the resistivity p. Shown is p 
vs. temperature, at AJaf=0, A=l, and Jaf equal to 0.14 (squares), 
0.16 (circles), 0.18 (closed triangles), and 0.2 (open triangles). The 
first two Jaf cases correspond to the FM regime, and the other two 
to the CE state. A first-order transition separates the metal from the 
insulator at low temperatures, with concomitant abrupt changes in 
the resistance, (a) corresponds to zero magnetic field pH/t=0 while 
(b) is at pH /i=0.05. Note in the latter that an insulating line in (a), 
turned metallic in (b), producing a huge MR. 

The most important result in this context occurs when mag- 
netic fields are introduced. In this case, the insulating behavior 
found in the CE regime close to the FM phase, more precisely 
at Jaf=0. 18, turns metallic in a first-order transition upon the 
application of a field of value 0.05£ (much smaller in magni- 
tude than the natural units of the problem, such as t). Assum- 
ing t of the order of 1,000 K, this is a field of about 50 T, larger 
but not by a huge amount with respect to those typically used 
in CMR experiments. By selecting Jaf even closer to the 
first-order transition the value of the field needed to induce 
metallic behavior can be easily reduced. The associated MR 
ratios - defined as MR=(p(0)-p(H))/p(H)xl00%, with H 
the field used - are shown in Fig.^] For the coupling where 
the insulator-metal transition was generated by the magnetic 
field, the MR ratio was found to be as large as 10 7 % at low 
temperatures, in excellent agreement with the values that can 
be deduced from the experimental data in Fig. ^] As ex- 
plained before, this drastic effect is simply caused by a level 
crossing induced by the magnetic field, which favors the FM 
phase. Note that the Jaf coupling used is not abnormally 
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close to the original transition, namely there is no need to 
carefully fix parameters to see this effect. The present cal- 
culation clearly shows that huge MR effects can be obtained 
in theoretical calculations, even in small clusters and without 
the need of tuning couplings. 
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FIG. 16: Clean-limit magnetoresistance (as defined in text) vs. tem- 
perature considering pH/t =0.025 (a), and pH/t=0.05 (b). The 
same set of parameters and convention of symbols are followed as 
in Fig. 1151 Huge MR ratios are observed, as found in experiments. 

Regarding the theoretical understanding of the more stan- 
dard CMR (namely, CMR2 in the present discussion), pre- 
vious investigations have relied on sophisticated calculations 
using simplified models and resistor networks in order to cap- 
ture the complex percolative physics expected to dominate in 
real compounds^. These subtle effects cannot be studied in 
the small clusters currently accessible to nearly-exact Monte 
Carlo studies of realistic models, as those presented here. 
However, results in those small clusters when quenched disor- 
der is included, already provide hints of the physics found in 
experiments and in simulations of toy models. For example, 
in Fig. I17h it is shown that the disorder reduces the resistiv- 
ity in the CE regime, and that this effect is further magnified 
when magnetic fields are applied. Unfortunately, in spite of 
observing these reasonable effects, the shape of the resistiv- 
ity curves still do not show the well-known metal-insulator 
transition at a finite temperature. Further computational work 
will be needed to reproduce the CMR2 effect using realistic 
Hamiltonians. Fortunately, the essence of the phenomenon 
appears to have been captured by the calculations of Ref.|5l 
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FIG. 17: Study of the resistivity p in the presence of quenched disor- 
der. Shown is p vs. temperature, for pH/t=0 (a), and pH /i=0.025 
(b), using A=l, A Jap=0.04 , and Jap equal to (squares), 0.04 (cir- 
cles), 0.08 (closed triangles), and 0.12 (open triangles). Large MR 
effects are observed, even involving only insulating phases. 



IX. CHARGE-ORDERING IN ELECTRON-DOPED 
MATERIALS 

For completeness, here results at electronic densities (n) 
different from 0.5 are also presented. However, the 'half- 
doped' character of the investigation remains. To be more 
specific, the emphasis in this section is on (n)=1.5, namely in 
the regime of electron doping of undoped compounds such as 
LaMnC>3, opposite to the hole-doping regime of (n)=0.5 dis- 
cussed in the rest of the paper. The experimental motivation 
for this effort relies on recent investigations that have reported 
results for Lao.7Ceo.3Mn03, which is indeed an electron- 
doped compound having LaMnOa as the parent material^. 
A ferromagnetic metallic phase has been reported in this con- 
text with a Curie temperature ^250 K and a large magnetore- 
sistance, establishing clear similarities with the well-known 
results for the hole-doped region of the phase diagram. 

These experiments suggest that there must be a relation be- 
tween electron and hole doping, which has not been discussed 
theoretically before to our knowledge. In particular, it is in- 
teresting to speculate what kind of states will be obtained in 
the electron-doped regime once the doping is made as large 
as 50% (i.e. (n)=1.5). In other words, here the study will 
focus on the discussion of what kind of CO/OO state could 
be obtained in the electron half-doped regime of manganites, 
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justifying the inclusion of the results described below in the 
present paper, which is devoted to half-doped manganites in 
general. 

§ # 4* # 

FIG. 18: Orbital, spin, and charge arrangement for the electron over- 
doped «n) = 1.5) CE/CO/OO phase. The 3a; 2 - r 2 and 3y 2 - r 2 
orbitals represent the Mn 3+ ions. The superposition of the x 2 — y 2 
plus 3z 2 — r 2 orbitals (closed circles) represents the Mn 2+ ions. 

The main conclusion of the analysis of the (n)=1.5 phase 
diagram using Monte Carlo simulations is that in the large 
Hund coupling limit, there is a mapping between (n) bigger 
and smaller than one (with (ra)=l corresponding to LaMnO.3). 
In fact, in the Jh=oo limit studied here the phase diagrams 
in the plane A- Jaf are identical (within numerical accuracy) 
with regards to the location of the phase boundaries. However, 
the characteristics of each phase cannot be the same in view 
of the different electronic densities. As particular case, our 
effort predicts that the analog of the well-known CE phase of 
half hole-doped manganites is the state shown schematically 
in Fig.^] In this state, there is a separation of charge into a 
'3+/2+' like configuration, to be contrasted with the standard 
3+/4+ state at (n)=0.5. The zigzag chains of the standard CE- 
state are still present, and even the orbitals in the middle of 
the segments of the zigzag chains are elongated along the x 
and y axis. The main difference between (n)=0.5 and 1.5 is 
the state at the vertices of those chains, which is empty at den- 
sity 0.5 but it contains two electrons at 1.5. An analog situa- 
tion is found in the other charge-ordered phases. It would be 
very interesting to confirm experimentally whether this state is 
present in real half-electron-doped manganites. Note that our 
calculation does not include the effect of Coulombic interac- 
tions, which may be more important in the electron-doping 
regime than in the hole-doping region. However, in other 
materials such as the cuprates, a strong similarity between 
electron- and hole-doping has been found as well, suggesting 
that a similar situation may occur in manganites. Complet- 
ing the phase diagram of electron-doped manganites would 
help us in achieving a deeper understanding of these materi- 
als. The concrete prediction of our efforts is that there should 
be an approximate particle-hole symmetry with respect to the 
undoped LaMn03 compound, in a similar spirit as it occurs 
in the cuprates™. 



X. CONCLUSIONS 

The investigations reported in this paper have unveiled sev- 
eral unexpected properties of half-doped manganites. For ex- 
ample, the CE phase was found to be more sensitive than ex- 
pected to the addition of quenched disorder. The effort in this 
context included the calculation of the clean-limit phase dia- 
gram in the T- Jaf plane, to understand the competition be- 
tween the FM and CE states believed by most manganite ex- 
perts to be at the heart of the CMR phenomenon. In agree- 
ment with previous calculations 2 ^, a low-temperature first- 
order transition between the two phases was found, with a 
shape similar to those reported experimentally^ Adding dis- 
order, the predicted^ depletion of critical temperatures was ob- 
served, but contrary to those expectations the reduction was 
not symmetrical in magnitude for the two competing phases. 
In fact, it was observed that the CE state rapidly transforms 
into a "CE glass" with disorder, while the ferromagnetic phase 
is comparatively less affected, in nice agreement with recent 
experiments where disorder is introduced "by hand" upon a 
clean-limit bicritical-shaped phase diagramliA^. The zigzag 
chains of the CE state and its very subtle arrangement of spins, 
charge, and orbitals seems easily broken by imperfections, 
contrary to the more robust uniform FM state. 

The present results revealed other new properties of man- 
ganite states that also deserve further investigations. Among 
these interesting new properties are the existence of novel 
phases, such as a FM/CO state, and a discussion on charge 
disproportionation that suggests that the widely held view of 
charge-ordering as containing 3+ and 4+ ions needs consider- 
able revision. 

A second important result presented here is the existence of 
a very large MR effect when the insulating state is close in en- 
ergy to the FM state, as it occurs in the first-order transitions 
systematically found in the clean limit. As a consequence, 
without quenched disorder and working at low temperatures, 
the insulating state can be destabilized by the FM state with 
increasing fields. The magnitude of the field required can be 
very small in units of the hopping, and its actual critical value 
depends on how close the CE state investigated is to the FM 
state in the phase diagram (i.e. how close the analyzed value 
of Jaf is to the critical value separating the phases). This 
MR effect is qualitatively different from the more standard 
effect found above the Curie temperature, which is believed 
to need disorder effects to be understood 2 ^. As predicted 
by Tokura and collaborators 27 many years ago, there are two 
types of CMR : a low-temperature form - investigated theoreti- 
cally here - and a higher temperature variety with the standard 
profiles for resistivity vs. temperature. In carrying out these 
investigations, we have introduced techniques borrowed from 
mesoscopic physics, within the context of the Landauer for- 
malism. 

The study of electron-doped systems has also been ini- 
tiated in the investigations reported in this paper, focusing 
on the half-doped limit. Charge-ordered states analogous to 
those observed at half-hole-doping were identified. A formal 
particle-hole symmetry appears to exist in the system with re- 
spect to the undoped LaMn03 limit, at least for a large Hund 
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coupling. More theoretical work in this area should be pur- 
sued, since experiments are starting to investigate electron- 
doped manganites. 

It is clear that the study of models for manganites us- 
ing unbiased techniques provides a plethora of interesting re- 
sults. In the field of manganites the crossfertilization theory- 
experiments has been remarkably positive, and new surprises 
will likely be found in the near future. Investigations of man- 
ganites and other related oxides should continue at its present 
rapid rate, to enhance our understanding of correlated electron 
systems in general. 
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XII. APPENDIX: SIMULATIONS IN THE CANONICAL 
ENSEMBLE 

In most of our simulations, we have kept the number of 
particles fixed in our systems basically at each Monte Carlo 



Step per Site (MCS/S). In order to achieve that, we have to 
solve the equation, 



N 2N 

EE r 

i=i k=i 



c Li c k,i 



cxp[(£ fe - fi)/T] 



for the chemical potential, /j,; where n is the desired number of 
electrons in the system, N is the number of sites of the cluster, 
cfj and Ck,i are the fermionic operators after diagonalizing the 
hamiltonian matrix (Eq. (1)), and Ek are the electronic energy 
levels. These energies are ordered in the way: E^n < Ek < 
Ei, with k running between 1 and 2N. The first sum symbol 
runs over sites, while the second runs over energy levels (2N 
in our case, since we have two orbitals per site). 



We have solved this equation using the Newton-Raphson 
method 31 , starting with an initial seed /xo = l/2(E2N-i + 
E2N-1+1), with I =int(xN), and the filling x given by: 
x = n/N. The symbol int( ) means the nearest integer of 
a real number. Fixing an absolute error for n of 10~ 5 , the 
chemical potential fj, is typically found in 4-5 iterations. At 
sufficiently low temperatures, k-sT/t <0.01 for example, the 
convergence of this method can fail. In this case we use the 
auxiliary bisection method 31 , where typically 20-30 steps are 
needed for the same accuracy. 
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